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I review the current status of the application of effective field theory to nuclear physics, 
and its present implications for nuclear astrophysics. 



1. Introduction 

> 

OO . 

l/-) \ When we finally arrive at the complete effective field theory that faithfully reproduces 
QCD in the low-energy regime relevant for nuclear physics, we will be able to answer 
some fundamental questions. In addition, to being able to describe inelastic processes, 
such as 3a — > 12 C + 7, in the same framework as elastic processes and the framework 
with which the nuclear energy levels are computed to high accuracy, we will know how 
such processes depend upon the fundamental constants of nature, the strong interaction 
mass-scale, Aqcd, and the masses of the quarks, m q . A small but significant step has been 
taken in this direction by a recent calculation of the behavior of the two nucleon systems 
as a function of m q . There are three good reasons for understanding such behavior. The 
first reason is intellectual curiosity. If we believe that we completely understand the strong 
interaction we should be in a position to address this issue. The second reason is to put 
constraints on physics beyond the standard model. There are hints that the constants 
of the standard model, such as a em and m g , may be time-dependent Jl.J. The third and 
perhaps the most practical reason is that the m^-dependence of nuclear properties and 
reaction rates needs to be known in order to extrapolate lattice-QCD calculations from 
the unphysical values that will be used in present simulations and in those of the near 
future, down to those of nature. For the foreseeable future lattice QCD calculations will 
be performed with unphysically large quark masses simply due to the fact that the time 
required to perform the simulations increases as one reduces the m q . Therefore, effective 
field theory (EFT) will play a key role in any future lattice QCD program, as it will only 
be through matching onto the appropriate field theory that lattice QCD will be able to 
make predictions of physical observables. This is true, of course, only until such a time 
when such observables can be computed directly with the physical values of the quark 
masses. The link between EFT and data has been and continues to be very strong, with 
data constraining the finite number of counterterms that appear at any given order in 
the EFT expansion. However, EFT calculations are being performed at sufficiently high 
orders so that in some cases the number of counterterms is greater than the number of 
observables, and lattice-QCD will be the only way to further improve the calculation. 
At this time constraints are being placed on linear combinations of the Gasser-Leutwyler 
coefficients, the Li, through partially-quenched QCD (PQQCD) simulations 0. 
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For processes involving multiple nucleons it is convenient to describe different momen- 
tum regimes with different EFT's. If one is dealing with processes in which all momenta 
are less than the pion mass, m^, it proves to be useful to use the pionless EFT, EFT (ft), 
while for processes involving momenta greater than m n one must use a pionful theory. 



2. The p <C Regime 

In the kinematic regime where |p| <C m v we can construct an EFT to describe the 
interactions of multiple nucleons and photons quite simply |3l4l5j . The fact that there 
is a bound state near threshold in the 3 S\ — 3 D\ coupled-channels, and a pole on the 
second-sheet near threshold in the ^So channel means that at least one operator in the 
EFT Lagrange density must be treated non-perturbatively. One is free to choose which 
operators are resummed, and this will be defined in the power-counting associated with 
the EFT. There are no explicit 7r's, p's or other hadrons that can be produced in the 
low-energy collisions. As EFT(;f) is only applicable in the very-low energy regime, chiral 
symmetry is not a good symmetry, however isospin is a good symmetry. The only input 
into the construction of EFT (ft) is Lorentz invariance, electromagnetic gauge invariance, 
baryon number conservation and the approximate isospin symmetry, the breaking of which 
is included perturbatively. 

The strong interactions between two nucleons are described by a Lagrange density 



C = Cn 



N T P j N 



f N T P j N + ... , (1) 



where P J is a spin-isospin projector, and the ellipses denote operators involving more 
powers of the external energy. It is clear that there is a correspondence between the 
operators in Eq. ((TJ and the coefficients in the effective range (ER) expansion of the scat- 
tering amplitude. This Lagrange density is explicitly Galilean invariant and relativistic 
corrections can be included in perturbation theory straightforwardly Further, it is 
straightforward to include the interactions in higher angular momentum channels, such 
as in the 3 S\ — 3 D\ coupled channels j^j- If the only thing that had been accomplished 
was to recover effective range theory from EFT (ft), nothing would have been achieved. 
However, EFT (ft) is more than effective range theory. In EFT (ft) one can systematically 
incorporate gauge fields, such as the photon and electroweak fields. In addition to interac- 
tions arising from gauging derivatives in the strong interactions, there are operators that 
are gauge-invariant by themselves. For instance, the lowest order (LO) gauge invariant 
operator that contributes to the deuteron quadrupole moment is 



Co = -eC(« 



N T RN 



N T PjN 



vV — l —y 2 ^ )A , (2) 
n — 1 



where Aq is the time-component of the electromagnetic field. A computation of the 
deuteron quadrupole moment up to next-to-leading order (NLO) allows one to fit the 
coefficient C^'. Once this constant is determined other observables that receive contri- 
butions from this E2 operator can be computed up to NLO [5]. 

The cross section for np —>■ at energies relevant to big-bang nucleosynthesis provides 
an example of a high precision calculation in EFT (ft), as shown in Fig. ^ This is a 
classic nuclear physics process in which meson-exchange currents play a significant role. 
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From the EFT(^") point of view, where mesons are not an explicit degree of freedom, such 
contributions are included via local gauge-invariant four-nucleon-one-photon interactions 
(e.g. with coefficient L np for np — > d'j). In this case the resonance saturation hypothesis 
works well, in so much as this local operator is saturated by one pion exchange (OPE), 
but in general this need not be the case. 




E T (MeV) E N [MeV] 

Figure 1. The left panel shows the cross section for 'yd — > np. The curves correspond 
to L np determined by the cross section for cold np — > d'j jH] and the dashed lines denote 
the ~ 3% theoretical uncertainty. Rupak has further reduced this uncertainty to below 
The right panel shows the total cross section for na scattering as a function of the 
neutron kinetic energy The diamonds and black squares are data. The dashed and 

solid lines show the EFT results at LO and NLO, respectively. 

Another place where EFT(^f) has had impact is in the analysis of SNO data, and the 
determination of solar neutrino fluxes. An integral part of the SNO analysis is the cross 
section for neutrino induced deuteron break-up, such as the charged current process v e d —>■ 
e~pp. In addition to the well-known one-body contributions to the break-up, there is also 
a contribution from two-body interactions. Butler, Chen and Kong [7] showed that only 
one constant, L^a, is required to describe the two-body contribution in order to calculate 
the cross-section at the percent level. Therefore, the question is where does one get L^a 
from? So far there have been three independent determinations of L^a- from the /3-decay 
of tritium in a hybrid EFT, from an analysis of reactor data on Vd break-up, and from 
an analysis of SNO data [HE!. They find L M = 6.5 ±2.4 fm 3 , 3.6 ±5.5 fm 3 , 4.0 ±6.3 fm 3 , 
respectively, which are all consistent within errors. 

Recently, Bertulani, Hammer and van Kolck have examined the relatively simple 
na "Halo-nuclear" system. They explored the cross section for na — > na by treating 
the na interaction as a sum of local interactions between a fundamental n-field and a 
fundamental a-field. There are three partial waves that make a significant contributions 
in the low-energy regime, J w = + , 1~, 1 + . The LO calculation requires the scattering 
lengths in each channel, while the NLO calculation requires both the scattering length 
and effective range in each. The experimental values a + = 2.4641 ± 37 fm, r + = 
1.385 ± 41 fm, ai- = -13.821 ± 68 fm 3 , n- = -0.419 ± 16 fm -1 , a 1+ = -62.951 ± 3 fm 3 , 
and rj+ = —0.8819 ± 11 fm -1 produce the curves shown in the right panel of Fig. ^ 

There has been continued progress in the three-body sector, and I will discuss only a 
small fraction of it here. Bedaque, Griefihammer, Hammer and Rupak JT] have refor- 
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mulated the EFT construction in the triton channel and have produced results for the 
nd — > nd phase-shift, as shown in Fig. |21 At LO and NLO there is a momentum inde- 




k [MeV] k [MeVI 

Figure 2. The neutron-deuteron 2 Si phase shift at LO (dotted line), NLO (dashed line) 
and NNLO (dark band) [H]. Left: Real part. Right: Imaginary part. The band in the 
NNLO curve results from varying the cutoff from A = 200 MeV to A = 600 MeV. The 
dots are the available direct phase shift analysis of the triton channel ^2] and the crosses 
are the results obtained with the Argonne V18+Urbana IX two and three body forces |13j . 

pendent three-body interaction, with coefficient H , that can only be determined in the 
three-body sector (or higher-body sectors), and the scattering length in the triton channel 
is sufficient. At NNLO there is a contribution from a momentum dependent three-body 
interaction, with coefficient H2, and the triton binding energy is sufficient to determine 
it. One sees that the EFT calculation is converging to the data and at this order shows 
no serious deviation from a more traditional potential model calculation. One can also 
examine the behavior of the Phillips line order-by-order in the EFT expansion. If one 
takes an arbitrary model that describes to arbitrary precision the two-body sector and 
uses it to compute the scattering length in the J = \ channel and triton binding energy, 
one will generate values that lie (approximately) on a line in the plane formed by them. 
This is the Phillips line, and results from the fact that if three-body interactions are 
possible then, in general, they will be present. The Phillips line, as shown in Fig. El is 
generated by varying the coefficient of the momentum independent three-body force over 
all possible values. 



a. (fra) 




B3 (MeV) 



Figure 3. The Phillips line [TI]. The dots correspond to the predictions of different models 
with the same two body scattering lengths and effective ranges The dotted and full 
line are the EFT results at LO and NLO, respectively. The cross is nature. 
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There have also been some exciting results from Hammer regarding the hyper-triton 
in EFT(^) ^5]. As the deuteron has only a small binding energy and the hyper-triton 
has a very small binding energy, B\ t = 2.35 ± 0.05 MeV, this particular three-body 
system can be described very well in EFT(f), and Hammer obtains a^ d = 16.8tl A fm and 
r^d = 2.3 ± 0.3 fm. I refer the reader to Ref. ^5] for a detailed discussion. An improved 
experimental determination of B^t would be welcome. 

3. The p ^> Regime 

It was Weinbergs' pioneering efforts ^H] i n the early 1990 's that led to the field of EFT 
in nuclear physics. He was attempting to construct an EFT for nuclear processes and 
nuclei involving momentum all the way up to the chiral symmetry breaking scale A x , 
and necessarily included the pion as a dynamical degree of freedom. The power-counting 
that he developed, known as Weinberg power-counting (W), involved performing a chiral 
expansion of the nucleon-nucleon potential using the same power-counting rules that apply 
in the meson and single nucleon sectors. The chirally expanded potential is then inserted 
into the Schroedinger equation to determine observables, such as phase shifts. However, 
there is a formal problem with this power counting jTTj in some of the scattering channels, 
particularly the ^So channel. However, extensive phenomenological studies with W power- 
counting appear to be in good agreement with data |18|19j . where such comparisons are 
possible, and the formal problem appears to have little impact when a massive regulator 
is used with a mass scale that is not radically different from a few hundred MeV. This 
problem led to KSW power-counting JT: in which the momentum independent four-nucleon 
operator is promoted to one lower order in the chiral expansion, and consequently pion 
exchanges are subleading order and treated in perturbation theory. This power-counting 
is formally consistent and gives renormalization group invariant amplitudes order-by-order 
in the power-counting. However, Fleming, Mehen and Stewart j2Hl (FMS) showed that 
the scattering amplitude in the 3 S\ — 3 Di coupled channels diverges at NNLO at relatively 
small momenta and KSW power-counting fails. FMS found that a contribution that 
remains large in the chiral limit destroys the convergence: it is the chiral limit of the 
tensor force that "does the damage" . Recently, it was suggested that one should expand 
observables about the chiral limit [21] (BBSvK power-counting). BBSvK power-counting 
has all the nice features of W and KSW counting: the chiral limit of the tensor force is 
resummed at LO along with the momentum and m q independent four-nucleon operator 
in the 3 5*i — 3 Di coupled channels, while pions are perturbative in the ^So channel and in 
higher partial waves where analytic calculations are possible. 

As I mentioned previously, the phenomenology of W counting has been explored ex- 
tensively. Epelbaum et al [22] h ave performed an impressive analysis of the few nucleon 
systems with W power-counting and find the calculations in the NN sector to converge 
well. The S-wave phase shifts are shown in Fig. HJ One sees that higher order calculations 
move closer to the data in all partial waves and the uncertainty is reduced. This calcu- 
lation uses a momentum space regulator, A, and one finds some sensitivity to the value 
chosen for A. An estimate of the uncertainty in this calculation can be made by varying 
A between 500 and 600 MeV, however, it would appear that such variations in A tend to 
underestimate the true uncertainty, as the NNLO result does not lie entirely within the 
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Figure 4. The S-wave phase-shifts for laboratory energies below 200 MeV |22| . Left 
(right) panel is the result at NLO (NNLO*). The momentum-space cut-off is chosen 
between 500 and 600 MeV leading to the band. The filled circles correspond to the 
Nijmegen PSA results [23] . 



NLO band. The convergence of the theory has also been examined in detail by Entem 
and Machleidt In other work Epelbaum et al j2H] have performed a detailed study of 
the three-body sector with W power-counting. One finds that the uncertainty is reduced 
by going from NLO to NNLO as expected, and their results are very encouraging. 




Figure 5. 4|'s for Compton scattering on deuterium [21]. The data are from Illinois [27] 
(circles), Lund [29J (diamonds) and SAL (squares) The solid line is the 0(Q 4 ) 

calculation with oln = 9.0 x 10 -4 fm 3 , (3n = 1.7 x 10 -4 fm 3 . The dashed line is the 
(parameter-free) 0(Q 3 ) calculation. 



Beane et al have obtained some exciting results for Compton scattering from the 
deuteron |2Sj- By working at order Q 4 in W power-counting the electric and mag- 
netic isoscalar polarizabilities of the nucleon were determined directly from the Compton 
scattering cross section at various energies, as shown in Fig. [2Sj. They find that 
a at = 9.0 x 10 -4 fm 3 , (3^ = 1.7 x 10~ 4 fm 3 . A higher order calculation is required in order 
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to reduce the uncertainties in a at and (3n and to ensure that the EFT calculation of this 
process is converging as expected. Further, it is clear that more data should be taken for 
this process, filling out the differential cross sections plots shown in Fig. 
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Figure 6. a + vs — a~ [HO]- The light and dark shaded regions are from the experimental 
pionic-hydrogen width and shift, respectively, taken from Ref. [HI]. The dotted line en- 
compasses the constraints from 7r-N phase shift data [32 . The dot is LO ;\TT. The two 
parallel bands are from ltd scattering [3U] evaluated with the NLO wavefunction with a 
cutoff of 500 MeV (upper curve) and 600 MeV (lower curve). 

Recently, W power-counting has been used by Beane et al to examine ltd scattering in 
EFT [30j. The EFT allows one to understand nd scattering in the same framework as ttN 
scattering. Given the precise data from pionic hydrogen and from pionic deuterium one 
can extract both the isoscalar and isovector ttN scattering lengths, a + and a~ respectively, 
with high precision. A plot of the constraints in the a + — a~ plane is shown in Fig. El 
from which it is found that or = 0.0918 ± 0.0013 m" 1 and a + = -0.0034 ± 0.0007 m" 1 . 
These values are in very good agreement with the recent analysis of Ref. [33] ■ A higher 
order calculation is required in order to further reduce the uncertainty in these values. 

An alternate approach to EFT calculations has been advocated by Rho and collab- 
orators. They argue that performing a chiral expansion of the current operators (or 
whatever operator you are interested in) and using wavefunctions generated by the best 
modern potentials ("best" being defined as those with lowest x 2 m t ne two- and three- 
nucleon sectors) should be equivalent to the formal EFT expansion. While I do not 
presently see how to formally justify this approach it appears to numerically converge 
well, and its cut-off dependence is systematically reduced as one performs calculations 
to higher orders. With this method, Park et al jH] have computed the cross sections for 
p p — > d e + v e and p 3 He — > 4 He e + v e . The wavefunctions are taken from the AV18/UIX 
potential while the current operator is expanded to a given order in W counting. To the 
order to which they work only one-body and two-body operators contribute. The one 
unknown coefficient appearing in the two-body contribution, L^a, is determined from 
the rate for tritium /3-decay, which they compute in the same framework. While the 
coefficient of the two-body operator depends sensitively on the value of the cut-off, the 
complete two-body contribution to the rate is relatively insensitive to the cut-off. They 
find Sh e p(0) = (8.6 ± 1.3) x 10~ 20 keV — b. This is a very encouraging result, as the 
calculation of Shep(0) suffers from significant cancellations. This method may provide a 
bridge between the formal EFT constructions and the many-body methods developed in 
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nuclear physics. 



4. The mg-Dependence of the Two-Nucleon Sector 



In some recent papers by Beane and myself [HI] and also by Epelbaum, Glockle and 
MeiBner j3E] EFT was used to determine the m q - dependence of scattering in the two- 
nucleon sector. In the ^-channel KSW power-counting can be used to derive an analytic 
expression for the scattering length, 



7 + 



g\M 1 



N 



**ft 



ml log 


(-) 







m. 



•) 2 -(7-^ 



M N ml 

47T 



{ 1 -^) 2 D 2 , (3) 



where 7 is a LO constant and D 2 (fi) is a combination of coefficients of operators with a 
single insertion of m q , that is presently unknown. The best that one can do at this point in 
time is to use naive dimensional analysis (NDA) to estimate a range of reasonable values 
for D 2 , defined by a parameter 77 <C 1 [HI]. The results of NDA are shown in Fig. [TJ NDA 




[MeV] m, [MeV] 

Figure 7. The left (right) panel shows the scattering length in the ^So-channel ( 3 Si- 
channel) as a function of the pion mass. The light gray region corresponds to r\ = 1/5 
and the black region corresponds to 77 = 1/15. In the 3 Si-channel the parameter die is 
taken to be in the interval -2.61 GeV~ 2 < 3 lfl < -0.17 GeV~ 2 and 3 18 = -0.51 GeV~ 2 . 



suggests that the di- neutron remains unbound in the chiral limit, while a relatively small 
increase in m q could lead to a bound di-neutron. 

In the 3 Si — 3 Di coupled channels the situation is somewhat more complicated. AT NLO 
in BBSvK counting not only does OPE contribute, but also the chiral limit of two-pion 
exchange. As a consequence, there are additional counterterms in the single nucleon sector 
that contribute in this channel but do not contribute to the ^So channel, in particular dig 
and die associated with the pion-nucleon interaction, and I4 associated with f n . This is 
in addition to the D 2 {ii) contribution in the 3 Si channel. The allowed regions for di$ and 
di6 are given in Ref.|JS], and I4 is known. Fig. [7| shows the presently allowed values of 
the scattering length in the 3 Si channel where we again have used NDA to estimate the 
possible values for D 2 . It is clear that for the range of parameters considered the deuteron 
could be bound or unbound in the chiral limit, and at present one cannot make a more 
definitive statement. This last statement disagrees with the conclusion of Ref. [33] . 
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5. Conclusions 

There has been rapid progress in the application of effective field theory to nuclear 
physics. A consistent power-counting has been established, improved model independent 
calculations of NN — * NN, Nd — > Nd, vd —>■ vd, pp — > dev, p 3 He — > 4 He ev, jd — > 7c? 
and red — > nd have been completed, and first efforts to describe halo-nuclei and hypernu- 
clei have been undertaken. Further, efforts to understand fundamental questions about 
nuclear observables have been made with interesting results. 

There are still many goals to be attained. It is clear that the calculations I have 
discussed here should be pushed to one higher order to convince ourselves about their 
precision. Moreover, it is clear that parameters in the single nucleon sector need to 
be known to higher precision than they are now in order to make precise statements 
in the multi-nucleon sectors. Ultimately, a much closer link between the lattice QCD 
community and the EFT community must be developed, as it is only with EFT that 
lattice calculations of complex hadronic systems will be possible, and further, it is only 
with lattice QCD that EFT will be able to make faithful calculations in multi-nucleon 
systems. This is clearly a symbiotic relationship that is yet to be realized. 



I would like to thank the organizers of PaNIC02 for putting together a stimulating 
meeting and inviting me to share the recent developments in EFT. I would also like to 
thank those who have let me reproduce their work in these proceedings. 
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